Nonuniform states in noncentrosymmetric superconductors 
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In noncentrosymmetric crystals, nonuniform superconducting states are possible even in the ab- 
sence of any external magnetic field. The origin of these states can be traced to the Lifshitz invariants 
in the free energy, which are linear in spatial gradients. We show how various types of the Lifshitz 
invariants in noncentrosymmetric superconductors can be derived from microscopic theory. 
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The discovery of superconductivity in CePtaSi (Ref. |l|) has renewed interest, both experimental and theoretical, 
in the properties of superconductors without inversion symmetry. One of the most spectacular differences from the 
usual, i.e. centrosymmetric, case is the presence of additional terms in the Ginzburg-Landau free energy, which are 
linear in spatial gradients. 2 These terms lead to a "helical" superconducting phase, in which the order parameter is 
nonuniform in the presence of an external magnetic field which is coupled only to the spins of electrons ££& That a 
nonuniform superconducting state can be created by purely paramagnetic effects was suggested a long time ago by 
Larkin and Ovchinnikov 6 and Fuldc and Ferrcll 7 (LOFF). In contrast to the helical state, the LOFF state appears as 
O , a result of the sign change of the second-order gradient term in the free energy. 

5^ ■ According to Refs. I3J3, the crystal symmetry might sometimes admit the linear in gradients terms - the Lifshitz 

Oh| invariants - in the free energy, leading to nonuniform superconducting states even in the absence of magnetic field. We 
^ ■ do not intend to provide a comprehensive list of possibilities for all noncentrosymmetric crystal symmetries. Instead, 
our purpose is to use two examples to show how the Lifshitz invariants at zero field can be obtained microscopically in 
noncentrosymmetric superconductors with a Rashba-type spin-orbit (SO) coupling. In Sec. U the case of a spin-triplet 
order parameter which transforms according to a three-dimensional irreducible representation of a cubic point group 
is discussed, while in Sec. |TT] we consider a mixture of two representations in a tetragonal crystal. Throughout the 
paper we use the units in which H = ks = lj and e denotes the absolute value of the electron charge. 
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, C -*, 1 I. INTERBAND PAIRING: SINGLE REPRESENTATION 

Our starting point is the following Hamiltonian of noninteracting electrons in a noncentrosymmetric crystal: 

(Si! Hq = J2 Yl M fc )<W + j(k)a a0 ]al a a k (1) 

O ■ k "0=U 
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where cr are the Pauli matrices, and the sum over k is restricted to the first Brillouin zone. The second term in Eq. 
(fT}, with 7(fc) = — 7(— fc), describes the Rashba-type (or antisymmetric) SO coupling of electrons with the crystal 
qq ' lattice, which is specific to noncentrosymmetric systems. In addition, there is a usual (symmetric) SO coupling, which 
is present even in centrosymmetric crystals. If the latter is included then a, j3 in Eq. (fT]) should be interpreted as 
pseudospin projections. The Hamiltonian is diagonalized by a unitary transformation a^a = X)a u a\(k)ck\, where 



- * — 1 

X 



UTA(fe) = Jhl+pi, UiX (k) = x —J^Jy — , (2) 



with the following result: 



ffo =EE^( fc )4AC fe A. (3) 



k X=± 



Here the band dispersion functions are £\(k) = eo(fe) + A|7(fc)|. The normal-state electron Green's functions can be 
written as 

G(k,w n ) = Yl tix(k)Gx(k,w n ), (4) 

A=± 



where 



Mk) = '-^^ (5) 



are the band projection operators (7 = 7/I7I), and 

G x (k,u>) = - ^TT (6) 

are the Green's functions in the band representation. 

In this section we show how one can obtain the Lifshitz invariants for an order parameter corresponding to a 
single irreducible representation of the crystal point group. Let us consider a purely triplet order parameter, which 
transforms according to an irreducible representation T (of dimensionality dp) : 

dr 

d(fc,r) = £>„(»•)<*,(*)> (7) 

a=l 

where tp a {k) = —cp a (—k) are the spin-vector basis functions, see Ref. llCl The order parameter matrix in the spin (or 
pseudospin) representation has the form A a p(k,r) = d(k,r)g a p : where g = ia&2- 

Using the standard formalism, see e.g. Ref. |4|, we obtain for the quadratic terms in the free energy density: 

F 2 = y E VaVa - Yl V* a K ab r] b , (8) 

a a,b 



where V > is the coupling constant, and the operator K a j, is obtained from 

d 3 k 



1 f d k 

Kab(q) = 2 T Ey (CT<i( fe Kj( fe )tr [#G(fc + q, u n )g G T (-k, -u n j\ 



(9) 



by replacing q — > D = — i V + (2e/c)A. The Matsubara summation is restricted by the cutoff energy u) c . Substituting 
here expressions ((4]) and using the identity 

tr [glfi Xl (k + q)g 3 fll(-k)} = 1 ~^ M 5 ij + A 1 A 27 »(fc) 7j (fc) - \{Xx - A 2 )e y - /7; (fe) 

(here we neglected the corrections of the order of q/kp), we obtain: 

1 r Mu 

Kab(q) = 2 T EEy ^<i( k )wAk){UWi(k)Gx(k + q, Un )G x (-k,-u> n ) 

+ [<>ij -li{ k )l]{ k ) - i^iji A ji{k)]Gx{k + q,w„)G_A(-fe, -w n )} . (10) 

The first (second) line of this expression describes the intraband (interband) pairing. 

The Lifshitz invariants originate from the odd in q contribution to the kernel, which in turn comes from the last 
(linear in 7) term in the second line of Eq. (p~0|) : 

d 3 k 

» A - <*0 

Neglecting the difference (of the order of |7|/ef) between the Fermi velocities in the two bands, i.e. setting v + (k) — 
V-(k) = vp(k), and calculating the Matsubara sums, we obtain in the vicinity of the critical temperature: 

Kab(q) = -i-^rq m eiji <$(fc)^ il (fe)^ J (fc) 7 i(fc)^, m (fc))^ , (12) 

where the angular brackets denote the Fermi-surface averaging, 



i I' d 3 k 

K^b(q) = -2^^ T EE A 7 j^<i^)<^Akmk)vx^(k)Gi{k,w n )G-x(-k,-u n ). (11) 



$(fc ) = Im^fl + ,MMU_^(3)|7( fc )l 
V 2 2t:T c o I TrT c o 



and ^(x) is the digamma function. We assume that the Rashba SO coupling is sufficiently weak in order for the 
interband pairing to survive. We note that in the absence of time-reversal symmetry breaking in the normal state, 
the basis functions can be chosen real. Then it follows from Eq. (fT2]l that the Lifshitz invariants are absent for order 
parameters transforming according to one-dimensional representations of the point group. 

Let us consider as an example a three-dimensional order parameter r) = (771, 772, %) in a cubic superconductor, which 
corresponds to the representation F\ of the point group G = O. We assume a spherical Fermi surface and describe 
the SO coupling by 

l(k) = 7 ofe, (13) 

where 70 is a constant. For the band dispersions we have £\(k) = eo(&) + Aa, where a — |7o|fc_F is the measure of the 
SO band splitting (kp is the Fermi wave vector). The normalized spin- vector basis functions have the following form: 

<Pa,i( k ) = y 2 e aijh- ( 14 ) 

Inserting this into Eq. (fT^j) . we obtain: K^ b (q) — i(7C(3)/8n 2 )(NFVFa/T^ )e a biqi- From this it follows that the 
Lifshitz invariant in the free energy density has the following form: 

F L = iK (r]lD y i] 3 + r)lD z rii + rj^D x rj 2 - c.c.) , (15) 

where 

8**!* v f { ' 

Note that, according to Eq. (fH)) . the order parameter Q satisfies d(k, r) _L ~f(k). In the band representation, this 
corresponds to interband pairing, as opposed to the limit of strong SO band splitting, in which only the component 
d || 7 survives (intraband or "protected" component), in addition to the spin-singlet component^ 

II. INTERBAND PAIRING: MIXTURE OF TWO REPRESENTATIONS 

In this section we discuss a different mechanism of producing the Lifshitz invariants. We start with the pairing 
interaction in the band representation, which can be written a general form as follows: 

F ™*=2v5Z Yl i A 2 (fe)iA3( fc ')^A 1 A 2 A3A4(fc I fe') c l+q,A 1 C -fc,A 2 C -fc',A3 C fe'+q,A 4 , (17) 

kk'q Ai.2,3,4 

where 

tx(k) = A 7x(fc)-»7,(fc) (18) 

are phase factors, 

Vx 1 A 2 A 3 A 4 (fc,fc') = V g (k, k')S\ 1 X 2 S\ 3 \ 4 +V u ,ij(k,k')Ti t \ 1 \ 2 (k)Tj t \ 3 \ i (k') 

+v m ,i(k,k')T i ,\ 1 \ 2 (k)6\ 3 \ i + v m ,i(k' ,k)8x 1 x 3 n,x 3 x i (k'), (19) 

and fj(fe) = u' (k)aiu(k), with the matrices u(k) defined by Eq. p]). The functions v g , v u> ij, and v m .i describe the 
pairing strength and anisotropy in spin-singlet, spin-triplet, and mixed channels, respectively. We follow the notations 
of Ref. [l2| and assume that the frequency dependence of the pairing amplitudes is factorized. The terms with Ai = A2 
and A3 = A4 describe intraband pairing and the scattering of the Cooper pairs from one band to the other, while the 
remaining terms describe pairing of electrons from different bands. 

Treating the interaction \Yl\ in the mean-field approximation, one introduces the gap functions A\ 1 x 2 {k,q) = 
t\ (k)A\ 1 \ 2 (k,q), with the following symmetry properties: 

AA 1 A 2 (fc,g) = A 1 A 2 A A2 A 1 (-fe,9). (20) 



These are related to the symmetry of the pairing interaction: V\ 2 \ 1 \ 3 \ 4 (—k,k') = \i\2V Xl \ 2 \ 3 \ A (k,k'), which can be 
easily established from Eq. (fTT|) . Near the critical temperature, the gap functions satisfy the linearized gap equations: 

/d 3 k' 
7^3 2J ^ a 2 a 4 a 3 (k, k')G X3 (fc' + q, w n ) G A4 (-fc', -w») Aa 3 a 4 (k',q), (21) 

n ^ ' A3A4 

where G x (k,uj n ) are the Green's functions of band electrons, see Eq. ((6|. 

We describe pairing anisotropy by the following model, which is compatible with all symmetry requirements: 

v g (k,k') = -V g , 

v Utij {k,k') = -V u (j(k)j(k'))S tJ , (22) 

Vm,i\ki k ) = U, 

where V g and V u are constants. Then, from Eq. (fTO|) one obtains: 

V XlX2X3Xi (k,k') = -V g 8 XlX J X3Xi - V u (j(k)i(k'))(T XlX2 (k)T X3X Jk')). (23) 

Further steps essentially depend on the crystal symmetry, which determines the momentum dependence of the SO 
coupling. Let us consider a tetragonal superconductor with the point group G — Ci V , in which case one can write 

7 (fc)= 7 o(zxfc), (24) 

where 70 is a constant. We assume a cylindrical Fermi surface along the z-axis. For the band dispersions we then 
have £ X {k) = eo(fc) + Xcx, where a = |7o|fc_F- From Eq. (|23p we obtain the pairing interaction components as follows: 

V++++ = V-— = -V g - V U M 2 , V++— = V—++ = -V g + V U M 2 , 

V+-+- = !>-+-+ = -V U M + KAf 2 , V+—+ = T>_++_ = -V U M - V U M 2 , 

V+++- = F + _— = V-+++ = F___+ = iV u MN, (25) 

V++-+ = V+-++ = V-+— = V—+- = -iV u MN, 

where M — kk' and N = (k x k') z . Note that the interband components in the second line here contain terms that 
are even in both k and fc', as well as ones that are odd in fc and fc'. As we shall see, this gives rise to inhomogeneous 
superconducting states. 

The pairing interaction can be presented as an expansion over the irreducible representations T of the point group 
G. The tetragonal group C^ v has four one-dimensional representations: Ai, A2, B±, B2, and one two-dimensional 
representation E (the notations are the same as in Ref. Il3h . The simplest polynomial expressions for the normalized 
basis functions have the following form: 

ip Al (k) = kl + kl = 1, ip Bl (k) = \/2(£; 2 - fc 2 .), <pB 2 (k) = 2V2k x k y > V>E(k) = V2(k x ,k y ). (26) 

Using 

M = -v> E {k)<p E {k'), 

M 2 = ~Lp Al (k)ipAi(k') + - V >B 1 (k)ifiB 1 (k') + ~ip B2 (k)tpB 2 (k'), 

MN = -^ Bl (fc)^ S2 (fc') - -<p B2 (k)<p Bl (k'), 
we obtain from the gap equations (|21[) the following expressions for the intraband gap functions: 

Aaa = V&PAiW +ifc<p Bl {k) + rfa<p Ba {k), (27) 

and also for the interband gap functions: 

A+_ = V *L<p Al (fc) + V+lv Bl (fc) + V+1<PbM + V+-fE(k). (28) 

A |_ can be obtained from A^ using Eq. (|2l7|) . The expansion coefficients r^ A , (q) play the role of the order parameter 

components. 



In a general nonuniform case and at arbitrary values of the SO band splitting and the coupling constants V g and V u 
[Eqs. l[2"2])]. the gap equations (p?Tj) yield a set of coupled equations for the eleven components of the order parameter. 
In order to demonstrate the possibility of inhomogeneous solutions even in the absence of an external magnetic field, 
it is sufficient to retain only the interband components in the A\ and E channels. To turn off the intraband pairing, 
we assume that V g < (i.e. the isotropic channel is repulsive), and introduce the notations rj,^ = r\ and f]^__ = £. 
Then, using the symmetry properties (|20j) . one can write 

A+_(fc, q) = Ti(q)<p Al (*) + t(q)>PE(k), A_+(fc, q) = -r]{q)<p Al (fe) + S(q)<p E (k). (29) 

The linearized gap equations take the form 

V = S^ V + S2^a, £a = S*TI + S%t b , (30) 



where a, b = 1, 2, and 



S™(q) =V U TJ2 f 4^3^ (k)G+(k + q, U»)GL(-fc, -to n ), 
S2HQ) = V U TJ2 / 7^3^(fe)^,a(fe)G+(fe + q,Un)G-(-k,-U n ), 

Sii(q) = V u TY,J^<PE,a(k)<PEAk)G + (k + q,uj n )G-(-k,-b Jn ). 

At q = 0, one obtains Effi = 0, which means that the A\ and E channels are decoupled. The critical temperature 
of the phase transition into a uniform superconducting state is the same for r\ and £: 

rp rpBcs ?C(3) a /q-i\ 

T cQ = T c0 - a , (31) 

16tt z 1 c0 

where T^ cs — (2e c /n)ui c ex.p(—l/NpV u ) is the standard BCS expression for the critical temperature in the absence 
of SO band splitting (i.e. at a = 0), ((x) is the Riemann zeta-function, C ~ 0.577 is Euler's constant, and Np is the 
density of states at the Fermi level (we neglect the difference between the densities of states in the two bands). It 
is assumed that a < T c q, otherwise the interband pairing is suppressed by the same mechanism that suppresses the 
singlet pairing in centrosymmetric superconductors. Indeed, the SO band splitting enters the Green's functions and 
the gap equations in exactly the same way as the Zeeman field does in the centrosymmetric case. 

The actual critical temperature of the superconducting transition is higher than T c q. Keeping the lowest terms in 
the gradient expansion in Eqs. p0[) . we obtain: 



[a(T - T c0 ) + Kg 2 ] 7? + Kqiti = 0, 

r k 

KqiT] + a(T - T c0 )% + — {5 lj q 2 + 2qiq 3 ) 



6-0, (32) 



where a = Nf/T c q, K — 7((3)Nfv f /32n 2 T 2 , and K = 2^/2Ka/vF- The linear in q terms, which mix the two 
channels, correspond to the Lifshitz invariant in the Ginzburg-Landau free energy. In the coordinate representation, 
this invariant has the following form: 

F L = k\rf{DZ) + (£*D)n], (33) 

where D — — iV, or, in the presence of magnetic field, D — — iV + (2e/c)A. It is easy to see that this expression 
satisfies all symmetry requirements. In particular, it is invariant under time reversal operation, which in the band 
representation is expressed as A\ 1 \ 2 (k) — ► A^ A (fe). According to Eq. (|29| . the order parameter components 
transform under time reversal as follows: r\ —* —7]*, £ — > £*. 
Seeking the order parameter in the form 

V (r) = Vo e^, $(r) = Coe^(l,0), (34) 

we obtain the following expression for the critical temperature as a function of q: 

5 




fl (T-T c0 ) = - i i^ 2 + -i^/32(-^) +rf (35) 



This function has a maximum at finite q = qo, where 



v F 
with c\ ~ 1.15. The corresponding critical temperature is 



qo = c\ — , (36) 



T c = T c0 + c 2 -(— ) , (37) 



where c 2 ~ 1.62. Thus, the nonuniform superconducting phase described by Eq. (134|) has a higher critical temperature 
than the uniform state. 

It is instructive to interpret our results using the spin representation of the order parameter: 

A a/3 (fc,q) = ^2 u *\i{k)A\ 1 \ 2 (k,q)u l3 x 2 {- k ) = -[u(k)A(k, q)u^(k)(ia 2 )] a /3- (38) 

A1A2 

The interband elements, see Eq. (f2"9"| . are translated into the spin representation as follows: 

A a p(k,q) = {iera 2 ) a pd(k,q), (39) 

where 

d x = -i-q{q)tp Al (k)k x , d y = -ir)(q)<p Al (k)k v , d z = -£(q)(p E (k). (40) 

Therefore, the pairing symmetry in our model is purely spin-triplet, with d(k, q) _L f(k). 

The Cooper pairs in the nonuniform state are composed of electrons from different SO-split bands, and the order 
parameter is modulated with the wave vector qo ~ (&/T c o)£o (here £0 is the coherence length). This effect formally 
resembles the nonuniform mixed-parity state (NMP) in centrosymmetric superconductors and superfluids in the 
presence of magnetic field, which was discussed in Refs. I14|)15L The reason is that, as was mentioned above, the SO 
band splitting in the noncentrosymmetric case affects the interband pairing of electrons of opposite helicity in the 
same way as the Zeeman field in centrosymmetric superconductors affects the usual BCS pairing between electrons 
of opposite spin. However, the NMP state originates from the triplet interaction channel introduced in the model 
along with the singlet channel. In contrast, the inhomogeneous superconducting state in the noncentrosymmetric case 
arises from a purely triplet pairing interaction, see the second line in Eq. (122[) . which acquires both fe-even and fc-odd 
components in the band representation. 

III. CONCLUSIONS 

We come to the conclusion that nonuniform superconducting states can exist in noncentrosymmetric superconduc- 
tors even in the absence of external magnetic field. Experimentally, the zero-field nonuniform states discussed in 
this article can only be observed in the noncentrosymmetric compounds with the SO band splitting smaller than the 
superconducting critical temperature. To the best of the authors' knowledge, in all noncentrosymmetric compounds 
discovered to date the relation between the two energy scales is exactly the opposite: the SO band splitting exceeds 
all superconducting energy scales by orders of magnitude, completely suppressing the interband pairing, both uniform 
and nonuniform. 

We would like to note that the models considered in Sees. HI and HT1 are different from those discussed previously in 
the literature. For instance, if the order parameter transforms according to a two-dimensional representation of the 
point group Dq (Ref. [8|), or has two intraband components, rj+ and ?y_, of the same symmetry in a tetragonal crystal 
(Ref. 9), then the Lifshitz invariant has the following form: 

F L = iK(r)* + D z r)_ - if-D zV+ ), (41) 

where if is a real constant. Although such terms satisfy all symmetry requirements and are therefore possible on 
phenomenological grounds, they are absent in our microscopic models. Moreover, expressions of the form (|4ip should 
exist even for the large SO band splitting, in which case the interband pairing channels are suppressed, and the 
mechanisms discussed in Sees. HI and HTl do not work. 
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